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ABSTRACT 



The goal of this research was to experimentally demon- 
strate the feasibility of a cylindrical shock wave to serve 
as the major pressure recovery mechanism for either 
cylindrical gas dynamic lasers or chemical lasers. An 
analysis was conducted to determine the minimum distance 
between the nozzle exit cylinder and the shock wave location, 
r m ~r e , required for power extraction. A cylindrical nozzle 
array was designed using the axisymmetric method of charac- 
teristics for a Mach number of 4.5 with a flexibility of 
other Mach numbers. A chamber incorporating end wall boundary 
layer bleed and ejectors for back pressure control was 
designed and built. The apparatus was run up to a pressure 
of 18 atmospheres. Three flow regimes were identified. A 

value of (r -r )/(r -r ) of 0.3 was attained, where r is 
m 6 S 6 s 

the shock wave location. Flow in the end wall nozzles 
consisted of a Coanda jet with a subsonic flow region 
adjacent to the wall; flow was not altered by boundary layer 
bleed. Except for the end wall nozzles a cylindrical shock 
wave was established. 
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I. INTRODUCTION 



At present there are three types of lasers that employ 
supersonic flow: chemical, CO electrical, and gas dynamic. 

All three are capable of producing large amounts of lasing 
power, but are restricted to use in platforms that are not 
limited by weight and volume. Due primarily to the large 
diffusers necessary to restore the flow to ambient pressure, 
large volume and weight are required to handle a given flow. 

A cylindrical geometry gas dynamic laser has advantages 
not found in the other supersonic lasers mentioned: 

(1) diffuserless operation 

(2) elimination of nozzle bending stresses 

(3) elimination of boundary layer corner interactions. 

Such a geometry would result in a compact system achieving 
lasing by rapid expansion in the radial direction rather 
than linearly as is presently done. The supersonic radial 
flow is then slowed to subsonic and raised to near ambient 
pressures by traveling through a cylindrical shock wave 
sufficiently downstream of the nozzles to allow for a 
reasonably large lasing cavity. A cylindrical shock wave 
eliminates the need for diffusers, thereby reducing the 
weight and volume of the device considerably. 

Supersonic conditions are achieved from flow traveling 
through axially symmetric nozzle rings, thus eliminating 
nozzle bending stresses found in linear gas dynamic lasers. 

A schematic of a cylindrical GDL is shown in figure (1) . 
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FIGURE 1 SCHEMATIC OF A CYLINDRICAL GDL (FROM 



Figure (2) shows schematic diagrams of two types of 
linear GDL's. Boundary layer buildup on the four walls of 
the lasing cavity interact at the corners, creating increased 
possibilities of flow separation and non-homogeneous conditions 
in the lasing medium. By having a flat plate at each end of 
the nozzle stack, the cylindrical geometry would eliminate 
such interaction. Since the radial flow experiences 
increasing area as it moves outward, the boundary layer buildup 
on each plate is less than in the linear case. 

A cold flow test of a number of nozzle rings was conducted 
to establish a cylindrical shock wave and to determine whether 
such an apparatus is feasible for use in a cylindrical gas 
dynamic laser. 

A theoretical power output example is included in section 
II of this report. Section III contains the development of 
the nozzle feed, the nozzle ring geometry and the ejector and 
boundary layer bleed systems. Experimental results are 
described in section IV. 
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(b) Open Cycle 



FIGURE (2) SCHEMATICS OF LINEAR GDL ' S 
(FROM REF. 5) 
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II. CYLINDRICAL LASER GEOMETRY 



A. PREVIOUS GDL’S 

Linear gas dynamic lasers are currently in use and have 
successfully produced laser beams in the order of a megawatt 
of power. These devices are quite large and heavy, thus 
limiting their use in operating environments (such as an 
aircraft) where such parameters are limited. These limita- 
tions are due primarily to the large diffusers necessary to 
return the supersonic flow to ambient pressure. 

A cylindrical GDL could reduce the weight and volume 
necessary to handle a given flow by expansion to supersonic 
velocities in the radial direction and returning to ambient 
pressure by means of a cylindrical shock wave. Other advan- 
tages include elimination of bending stresses and reduced 
boundary layer interferences. 

Ortwerth [Ref. (1)] designed a cylindrical source flow 
device, but retained diffusers similar to those used in 
linear GDL's. The diffusers were aligned radially as shown 
in figure (3) . Such a design reduces the volume required of 
such a device by having an annular lasing cavity rather than 
a linear one. Retaining the diffusers does not allow for 
the substantial reductions in weight and volume that a 
cylindrical shock wave would facilitate. 

Sponsored by the Air Force Weapons Laboratory, the AVCO 
Everett Research Laboratory [Ref, (2)] conducted an extensive 
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FIGURE (3) SCHEMATIC OF A CYLINDRICAL GDL WITH DIFFUSERS 



theoretical program to establish the feasibility of 
cylindrical source flow for a cylindrical GDL. Included 
in the program was a cold flow test of one nozzle ring which 
established supersonic flow at the predicted mach number but 
failed to create a cylindrical shock wave. Results showed 
that homogeneity of the lasing medium can be made 
acceptable. Separation of the boundary layer interfered with 
the establishment of a cylindrical shock wave, and it was 
recommended that some type of boundary layer control be used. 

Figure (4) shows the assembly used in the AVCO Everett 
test . 

B. THEORETICAL POWER OUTPUT 

In a linear GDL the diffusers can be placed far enough 
downstream so that a large lasing cavity is assured. In 
a cylindrical GDL the width of the lasing cavity is deter- 
mined by the shock location. To extract all the power 
available from a given flow, the cylindrical shock must not 
occur too near the nozzle exit. Figure (5) shows the nozzle 
stack with the lasing mirrors. 

Assuming a Gaussian intensity distribution on each mirror, 
the intensity I is given by 

peak 6 



I (r) = I 



(1) 
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FIGURE (4) AVCO EVERETT TEST RIG (FROM REF. 2) 
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FIGURE (5) CYLINDRICAL NOZZLE STACK WITH LASING MIRRORS 



Satisfying the boundary conditions and re-defining r yields 



Kr) = I p e 



-[£r£=VT ) 2 

r e -r kilowatts 

1 



( 2 ) 



cm 



where r g and r are shown in figure (A-l) in Appendix A. 

Integrating over the mirror area gives the total power 



m 

P = 2 I / r expt^-^-VT ] 2 dr 

p r r -r 

e e 



(3) 



An alternate method of finding the power available from 
the flow is 



Power = (mass flow rate) x (specific energy) 



or 



P = 2 r e Lp Q a o (p/p Q ) ( a e /a Q ) M g x (E) kilowatts (4) 

where subscripts "e" and "o” refer to exit and stagnation 
conditions respectively. Equating equations (3) and (4) gives 

r m 

I / r exp[£l£-V2 ] 2 dr = r e Lp a (p/p ) (a e /a ) M q E 
r r -r 



where the only unknown is the outer radius of the mirror, 
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Performing the integral in equation (5) and solving for 



r 



m 



r 



m 




( 6 ) 



where c 



4r e L Po a o (a/a o )M e ^P/Pq) 






Using the values for the constants listed in Appendix A, 
equation (6) reduces to 



where L is the overall length of the nozzle stack. Figure 
(6) shows the minimum shock location for increasing nozzle 
stack length. 

For six nozzle rings, the outer radius of the mirrors 
is 1.77 inches. Using this value in equation (3), the 
power is found to be 32 kilowatts. Details of the theoretical 
power available calculations are shown in Appendix A. 




(7) 
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FIGURE (6) OUTER RADIUS OF LASING MIRROR VS. NOZZLE STACK LENGTH 



III. EXPERIMENTAL DEVELOPMENT 



A. DESCRIPTION OF EXPERIMENT 

Reference 2 is a final report of analytical and experi- 
mental programs conducted by the AVCO Everett Research 
Laboratory to determine the feasibility of using a cylindrical 
geometry for a gas dynamic laser. The Naval Weapons Center, 
China Lake, California, expressed an interest in the concept 
of a cylindrical GDL and requested that research be conducted 
to determine whether a cylindrical shock wave could be 
established by incorporating the recommendations of the AVCO 
Everett report. Funding of this project was provided by 
Naval Air Systems Command, Code 320. 

Before determining the actual shape of the nozzle rings 
(see figures (8) and (9)), it was first necessary to calculate 
the throat height, that is, the spacing between each nozzle 
ring in the stack. This spacing is determined by the mass 
flow that could be supplied to the nozzle stack, which is 
limited by the feed holes. See figure (14) . 

Once the throat height was determined, the exit to throat 
area ratio, and thus the exit Mach number, was found from the 
geometry of the nozzle rings. The axisymmetric method of 
characteristics found in Ref. 3 was then used to determine 
the contour of the nozzle walls. Once the equation describing 
the wall shape was calculated, the design was sent to the 
machine facility at NWC , China Lake, where the nozzle rings 
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were fabricated. Figures (8), (9), (10), (12) and (14) were 

provided by NWC, China Lake. 

To prevent separation of the boundary layer at the shock 
location, suction was used to bleed off the low energy flow 
in the lower portion of the boundary layer. Extensive 
analytical investigation of supersonic boundary layer bleed 
was conducted and is found in section III-D of this report. 

The loss in total pressure across a normal shock of the 
strength anticipated resulted in total pressure downstream 
of the shock (using a reservoir pressure of ten atmospheres) 
to be less than one atmosphere. Ejector nozzles were 
necessary to lower the back pressure seen by the flow exiting 
the nozzle rings and thus allow them to start. Design of 
the ejector apparatus is found in section III-E of this 
report. 

Figure (7) shows the assembled experimental apparatus 
with dimensions. The design allowed for testing of up to 
eight nozzles at the design Mach number. The assembly shown 
was attached to the free jet located in building 230 of the 
Naval Postgraduate School. Air was supplied to the free jet 
from a 2400 cubic foot tank capable of pressure up to three 
hundred pounds per square inch. The moisture was removed 
from the air before it entered the tank by means of a 
hydryer. A model 06231 Lammert Industries vacuum pump was 
used to evacuate a 200 cubic foot tank to provide suction 
for boundary layer bleed. 
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FIGURE (8) ASSEMBLED NOZZLE STACK (COURTESY NWC , CHINA LAKE) 
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FIGURE (9) DISASSEMBLED NOZZLE STACK (COURTESY NWC , CHINA LAKE) 






B. NOZZLE FEED SYSTEM 
1 . Throat Height 

It was found, as expected, that the throat height 
was the most critical dimension of the entire design. A 
relatively large throat height is desirable in that the 
ratio of the nozzle exit area to the throat area is kept 
small, thus keeping the exit Mach number reasonably low. 

This would have the advantages of ease in starting and 
smaller total pressure loss across the shock. 

A large throat height also requires large mass flow 

requirements, and the amount of air that can be fed to the 

nozzle rings is limited by the inlet holes in the center of 

each ring. Figures (10) through (14) show the dimensions 

and photographs of the nozzle rings. The area of the feed 

holes must be greater than the combined throat areas of the 

nozzle ring stack to prevent choking of the flow at the feed 

holes. Thus if many nozzle rings are to be tested, an 

extremely small throat height must be used. This would 

result in a number of complications: 

(i) An excessively large A exit / A throat / causing 

high supersonic or hypersonic flow at the nozzle 
exit and thereby creating a problem in starting 
and a large loss in total pressure across the shock. 

(ii) Extreme heat concentrations at the throat (not a 
problem in cold flow testing) . 

(iii) Difficulty in fabricating spacing washers of uniform 
thickness, causing unacceptable variations in the 
flow due to changes in the throat height. 
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FIGURE (10) SIDE VIEW PHOTOGRAPH OF SINGLE NOZZLE RING (COURTESY NWC , CHINA LAKE) 
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FIGURE (11) DIMENSIONS OF NOZZLE RING (SIDE VIEW) 
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FIGURE (12) END VIEW PHOTOGRAPH OF SINGLE NOZZLE RING (COURTESY NWC , CHINA LAKE) 







FIGURE (13) DIMENSIONS OF NOZZLE RING (END VIEW) 
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FIGURE (14) END VIEW OF NOZZLE STACK 



The throat height must satisfy the continuity 



equation : 



m in 



6 

= E 
n=l 



m 



n 



where n is the number of nozzle rings. This can also be 
expressed as 



(2p . U. A. ) 
'in in in 



= 6 p*U*A* 



(1) 



Rewriting equation (1) using the equation of state and 
solving for the throat height: 



h = 



_ 1 _ 

6tt 



.. _ , P in. , T* . h 

M in A in "f^~ T. ^ 
in 



Using the value of inlet area previously obtained and using 
the isentropic relations for pressure and temperature in 
equation (2) , the throat height in inches can be expressed 
as 



.092 M 

(1 + * 2M in )3 



(3) 



T The factor of two in equation (1) results from the 
nozzle stack being fed from both ends. 
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Equation (3) shows that the throat height is a function of 
inlet Mach number only. Details of the throat height 
calculations are found in Appendix C. 

The exit Mach number is determined by the ratio of 
the exit area to the throat area. This ratio is found from 
the geometry of figure (11) to be 



A exit 

A 



2irr e (h+t) 
2Vr^h 



or 



A e _ 1 . 5 (h+, 5) 
A H 



(4) 



Also, assuming isentropic flow through the nozzle, the area 
ratio can be expressed as 



A 

e 

A 



_ 1 _ 

M 

e 







M e 2 )] 



Y+l 

2Ty-TV 3 



(5) 



Equations (4) and (5) show that the exit Mach number is 
directly dependent upon the throat height. Since the throat 
height is a function only of the inlet Mach number, the exit 
Mach number is also a function of the inlet Mach number. 

The results of equations (3) , (4) and (5) for various inlet 
Mach numbers are listed in Table (1) . Figures (15) through 
(17) show the results of Table (1) in graphical form. 
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Table 1 


Throat Height and Exit Mach 
Various Inlet Mach Numbers. 


Number for 


M inlet 


h (inches) 


M . . 
exit 


.2 


.0179 


5.65 


.4 


.0334 


4.95 


.6 


.0446 


4.62 


.8 


.0511 


4.47 
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FIGURE (15) THROAT HEIGHT VS. INLET MACH NUMBER 
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FIGURE (16) EXIT MACH NUMBER VS. INLET MACH NUMBER 
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FIGURE (17) EXIT MACH NUMBER VS. THROAT HEIGHT 



A throat height of .05 inches was chosen to allow 
for the lowest exit Mach number that would not result in a 
choked inlet condition. From figure (17), a throat height 
of .05 inches results in an exit Mach number of 4.5. The 
axisymmetric method of characteristics was used to design 
the contour of the nozzle rings. These calculations are 
found in section III-C. 

Knowing the throat height between each nozzle ring, 
the mass flow rate can be determined. For each nozzle, 



m nozzle 



p*U*A* 



Using the equation of state and rearranging yields 



P 0 (]H (49.1) 

* r O t 

m = [ ] 

T* h 

R(T 0 £-> 2 



. 2tt (1) (105) 

[ ] . 

144 



Substituing the appropriate values, the mass flow per 
nozzle is found to be 



= 1.08 lb m/sec 



Total mass flow is 



“total “ 6 (m n> 



m, . , = 6.48 lb /sec 

totctl m 
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2. Axial Variation of Mach Number and Static Pressure 



After the flow enters the nozzle stack, changes in 
Mach number and static pressure occur in the axial direction. 
Shapiro [Ref. 4] develops influence coefficients that relate 
changes in one flow condition to changes in another. 

Assuming constant molecular weight, the change in Mach 
number is related to the area change by 



2 

d.M z 




2(1 + I r- m 2 ) 

, „2 A ' 

1-M 



(6) 



Appendix D shows the model used for the axial flow through 
the nozzle stack. The pressure change is also related to 
changes in area and is given by 



dP _ yM 2 dA 

P ~ T A 

1-M 



Using the geometry of the nozzle stack, and rearranging 
equations (6) and (7) , axial Mach number and pressure 
changes are found to be 



1 + 1-1 M. 2 

M.., = M. - M. [ — — * — — — (.06667)] 

l+l i i i _m. 



( 8 ) 



and 

2 

P.., = P. + P. [ =-*■ (.06667)] . (9) 

l+l i i i-M, 
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A computer program was written using equations (8) and (9) 
to determine values of Mach number and pressure in the 
axial direction. Figures (18) and (19) are the results of 
this program in graphical form. A listing of the program 
and the results are found in Appendix D, along with details 
of the preceding development. 

3 . Radial Variation of Pressure and Mach Number 
As the flow leaves the nozzles, it experiences 
increasing area as it moves out radially. This area 
increase results in increasing Mach number and a corresponding 
decrease in static pressure. 

From the nozzle ring geometry, 



a/a* = 2 r r(h+t) 

A/A 2iTr, h 



A/A* = 11.0(r) 



(10) 



Assuming isentropic flow: 




(ID 



and 



Y 




( 12 ) 
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FIGURE (18) MACH NUMBER VS. AXIAL STATION 
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FIGURE (19) PRESSURE VS. AXIAL STATION 



Equations (10) , (11) and (12) relate Mach number and 

pressure to radial position. Table (2) is a summary of 
these three equations for different radii. Figures (20) 
and (21) are graphical representations of Mach number and 
pressure as functions of radius. 



C. NOZZLE RING DESIGN 
1 . Nozzle Contour 

As was shown previously in Section III-A, an exit 
Mach number of 4.5 was chosen for design. The curvature of 
the nozzle rings was determined by the use of the axisym- 
metric method of characteristics, 
a. Initial Turning 

In a two-dimensional nozzle the flow is turned 
through ( v fj_ na i/2) degrees at the throat, where v f^ na ^ is 
the value of the Prandtl-Meyer function corresponding to the 
design exit Mach number. In a three-dimensional case the 
flow expands from throat to tip (and therefore the Mach 
number increases) due solely to the increase in area. Thus 
to achieve the desired exit Mach number in the three dimen- 
sional case, the flow should be turned less than in the 
two-dimensional case. 

Assuming parallel plates, the ratio of the exit 
area to the throat area reduces to 



A e/ A throat “ 1 - 5 / 1 • 



Table 2 


Radial Variation of Mach 


Number and Pressure 


r (in.) 


14 


• p stat ( P sf) 


1.5 


4.50 


73.11 


2.0 


4.84 


48.30 


2.5 


5.089 


36.07 


3.0 


5.239 


30.39 


3.5 


5.479 


23.27 


4.0 


5.674 


18.85 


4.5 


5.869 


14.92 
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FIGURE (20) MACH NUMBER VS. RADIUS 



STATIC PRESSURE 
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FIGURE (21) PRESSURE VS. RADIUS 



Using a design Mach number corresponding to an area ratio 
of 



A/A* = (A/A* ) M = 4. 5 (1/1. 5) 



or 



A/A* =10.4 , 

and then treating the initial turning as if it were a 
two-dimensional case should give approximate results (that 
is, an exit Mach number near 4.5). For 



A/A* =10.4 , 



M . . 
exit 



3.96 



which yields 

v 

exit 



65.3° 



Therefore the initial turning of the flow at the throat 
should be 



A9 



v 

exit 
“1 — 



32.65° 



where 6 is measured from the flow direction entering the 
throat . 
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b. Axisymmetric Method of Characteristics 

Liepmann and Roshko [Ref. 3] develop expressions 
for finding the Prandtl-Meyer function (v) and the flow for 
direction (0) for an axially symmetric flow and are given by 



and 



v 3 “ 7 (v l + v 2 } + 7 (0 1 " 0 2 } 



^ sin 6. 

+ j[sin Ul — 



sin 0. 



AC 13 + sin v 2 _ 



An 23 ) 



(13) 



e 3 - r (v l - u 2> + l< e i + e 2> 



, sin 0. 

+ 7 [sin m ± — 



sin 0, 



A5 13 - sin y 2 — 



in 23 1 



(14) 



where all angles and subscripts are as shown in figure (22) . 

Using equations (13) and (14) the nonsimple core 
of the characteristic mesh was calculated. Figure (23) shows 
the entire characteristic mesh. The initial turning at the 
sharp-edged throat was divided into eight segments as figure 
(23) illustrates. Section II of Appendix E contains tables 
of Mach number, Mach angle, flow angle and Prandtl-Meyer 
function for each point in the mesh. 

As the last of the right running characteristics 
from the throat reached the plane of symmetry (point 70) , 
the Mach number was found to be 4.47, and the flow angle was 
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FIGURE (22) CHARACTERISTIC ELEMENT 



53 




54 



FIGURE (23) COMPLETE CHARACTERISTIC NETWORK 



90°. These flow conditions show that by reducing the 
design Mach number by the method previously described, and 
then treating the initial turning as a two-dimensional case, 
accurate results will be obtained. 

c. Establishing Points on the Nozzle Wall 

The left running characteristic from point 70 
is a straight line, and the flow conditions at any point on 
this characteristic will be the same as at point 70, that is, 
M = 4.47 and 9 = 90°. Points on the nozzle wall were deter- 
mined by calculating points along a right running character- 
istics in the -£ direction starting from a point on the left 
running characteristic between points 70 and 76. 

Following a development similar to Ref. 3, 
equations (13) and (14) were modified to calculate points 
of the characteristic mesh proceeding in the -E, direction. 

The modified equations are 



and 



v 3 = r (v l + v 2 ^ + j( 0 i “ e 2^ 



+ j[-sin \i 1 



sin 0 . 



AC 13 + sin y 



sin 0 ~ 

2 — — ai, 23> 



(15) 



0 3 = 7 (v l " v 2 } + 7 (0 1 + 0 2 } 



1 sin 0 

+ j-t-sin Ul - — ■ 



sin 0 , 



AC 13 “ sin y 



2 r. 



An 23 3 



(16) 
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where angles and subscripts are as shown in figure (E-l) 
found in Appendix E. Equations (15) and (16) were used to 
construct right running characteristics in the -£ direction 
from points 72, 74, 76 and 78 (not shown) until the left 
running characteristics reflected from the centerline were 
met. Starting at the throat the wall was drawn at 32.65° 
(relative to the centerline) until point 19 was reached, 
where the wall angle was reduced to the flow angle that was 
calculated for point 19 . This procedure was continued for 
points 32, 37 and 50. The last point (point 80) was the 
result of changing the angle of the wall to the estimated 
flow angle at that point. This line was extended until a 
radius of 1.5 inches was reached. Limiting the exit radius 
to this value was an effort to keep the exit Mach number 
as low as possible as was discussed previously. 

d. Wall Shape 

The r and z coordinates of the points where the 
wall slope changes, as well as the throat and the exit points, 
were fit with a fourth-degree polynomial using the Hewlett 
Packard model 9830-A computer and a polynomial regression 
program. The coordinates of the points and the coefficients 
of the fourth-degree polynomial are shown in Table (3) . 

Figure (24) shows the nozzle contour from throat to exit. 

2 . Trailing Edge 

As can be seen from figure (24) , the nozzle contour 
is not tangent to the ring centerline at r = 1.5 inches. 
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Table 3 


Description of Wall Contour 


A. Wall Point 


Coordinates from Method of Characteristics 



Point Number 


X Coord. 


Y Coord. 


(throat) 


1.00 


.025 


19 


1.06 


.066 


32 


1.11 


.097 


37 


1.18 


.138 


50 


1.29 


.173 


80 


1.44 


.215 


(exit) 


1.50 


.226 


B. Coefficients of 


Fourth Degree 


Polynomial Fit to Points 



B (0) = -6.375 
B (1) = 18.144 
B (2) =-19.773 
B (3 ) = 9.939 

B (4 ) = -1.910 
R-Square = 0.99997 
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FIGURE (24) NOZZLE WALL CONTOUR: FOURTH DEGREE POLYNOMIAL FIT TO CALCULATED 

WALL POINTS 



Rather than change the shape of the wall to a contour that 
would result in a sharp trailing edge, it was decided to 
maintain the design contour. This resulted in the flow 
exiting the nozzle at a 5° angle and a blunt trailing edge 
as shown in figure (25) . 

3 . Throat Shape 

A sharp-edged throat was used to obtain the highest 
degree of turning for the shortest length. Reference 5 
states that to be assured of uniform flow at the nozzle 
entrance, the radius of curvature of the inlet should be 
two to three times the throat height. Figure (26) shows the 
details of the nozzle inlet and throat. Due to the original 
geometry of the nozzle ring (before machining) the ratio of 
the radius of curvature to the throat height was limited to 
between 1.5 and 2.0. The shaded area in figure (26) 
represents the machining tolerance of the inlet contour. 

The center of curvature was displaced to the left of the 
sharp corner to assure horizontal tangency of the flow. 

D. EJECTOR NOZZLE DESIGN 

For a reservoir stagnation pressure of ten atmospheres 
and Mach-five flow, the total pressure downstream of the 
normal shock will be less than one atmosphere. To lower 
the back pressure and thus allow the cylindrical nozzles to 
start, ejector nozzles were designed and installed in the 
exhaust section of the apparatus as illustrated in figure 
(7). Using the continuity, momentum and energy equations. 
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FIGURE (25) DETAIL OF NOZZLE RING TRAILING EDGE 





FIGURE (26) DETAIL OF NOZZLE RING THROAT 
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solutions to the flow were found for different combinations 
of exhaust area and ejector nozzle diameter. A computer 
program was written to solve the three equations, and the 
results are tabulated in Appendix F. Figure (27) is a graph 
of these results. Any combination of nozzle diameter (Dl) 
and exhaust, width (D3) which results in a chamber pressure 
(PI) below the total pressure line will allow the cylindri- 
cal nozzles to start. A nozzle diameter of 1.75 inches and 
exhaust width of 6.5 inches were chosen. This combination 
results in a back pressure of .562 atmospheres. Complete 
details of this development are found in Appendix F. 

E. BOUNDARY LAYER BLEED 

Boundary layer bleed was employed to keep the flow from 
separating at the shock location. By bleeding off the low 
energy air in the boundary layer the remaining flow will 
have sufficient energy to support the pressure rise associ- 
ated with a shock of the strength anticipated. A search of 
the literature relating to boundary layer failed to locate 
specific information on supersonic boundary layer bleed 
techniques, specifically, hole size and pressure ratio 
required to remove a certain fraction of the boundary layer. 
Since information could not be found, a model was developed 
to make such a prediction. 

1 . Flow Model 

A model fluid- discussed by Stewartson [Ref. 7] was 
used in this development. The fluid is assumed to have a 
Prandtl number of unity and the viscosity is proportional to 
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FIGURE (27) BACK PRESSURE VS. EXHAUST WIDTH WITH EJECTOR NOZZLE DIAMETER 
AS A PARAMETER 



the static temperature so that up = U-^P^/ where subscript 1 
denotes the reference condition. Assuming the fluid has a 
Blasius profile and that stagnation temperature is constant 
through the boundary layer, the incompressible profile was 
transformed to a compressible profile by the method described 
in Ref. 7. 

For analytical ease, a Blasius profile was 
approximated by 

u/U = sin (An) 

where A = 54° and n is given by Ref. 7 as 

n = X V2/R e 

The reference length is denoted by X. The model shown in 
figure (28) was used to describe the flow. By applying a 
bleed plenum pressure (P2) lower than the freestream 
pressure (P^) , a streamline (4) is assumed to stagnate on 
the downstream edge of a slot of width H. All flow below 
this streamline goes through the slot, which is assumed 
long enough to have parallel streamlines at the slot exit. 

By assuming isentropic flow in each streamtube below (4) , 
it is possible to calculate flow conditions at the slot 
exit by specifying (P^/P 2 ) . Using the continuity equation 
leads to a slot width, H. The complete set of equations 
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programmed for the Hewlett Packard model 9830 computer, 
and is listed in Appendix G. 

2 . Results of Calculation 

Figure (29) illustrates the results of the computer 
program and plots the fraction of the boundary layer that 
enters the slot (Y/6) as a function of the pressure ratio 
(Pl/P 2 ) , with nondimensional slot width H/6 as a parameter. 
Table (4) gives the results for the flow conditions antici- 
pated for radial flow at a Mach number of 4.5. The symbols 
used are shown in figure (28) or have been given previously 

with the exception of a , which is the stagnation speed 

e 

of sound at the edge of the boundary layer. 

The dot-dash-dot line in figure (29) represents the 
locus of maximum Y/o for a given slot width and pressure 
ratio. 

3 . Bleed Hole Design 

Using the values from table (4) , (specifically H and 
X) the total bleed area was found to be 

A, , , = 2 ir (r + X) H 

bleed e 

= 2 ir (1 . 5 + 3) (.00143) 

^bleed = ,38 ln 

Three hundred holes were drilled through each plate at a 
diameter of .141 inches as shown in figures (30) and (31). 
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FIGURE (29) PLOT OF FRACTION BOUNDARY LAYER BLEED, Y/6, 

AS A FUNCTION OF PRESSURE RATIO AND NON- 
D I MEN S I ON AL BLEED SLOT WIDTH (FRuM REF. 8) 
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Table 4 Flow Conditions for Boundary Layer Bleed (from Ref. 8) 



M = 4.5 
y = 1.4 

U e = 2318 ft/sec 
6 = 0.00504 ft 
H = .00143 ft 

STREAMTUBE 

NUMBER n Y 

E-3 

ft 



0 0 



1 






.40 


.60 


.076 


1.56 


0.74 


0.15 




.16 


.40 














2 






.40 


1.88 


.227 


1.58 


2.21 


0.47 




.32 


.80 














3 






.39 


3.01 


.381 


1.63 


3.67 


0.77 




.48 


1.19 














4 






.38 


3.98 


.537 


1.69 


5.10 


1.04 




.64 


1.58 












- 


5 






.37 


4.75 


.698 


1.78 


6.51 


1.28 



.80 1.95 



R e = 500,000 



x = 0.25 ft 



A = 54° 

= 200 psf 
Y/6 = 0.385 



a Q = 1157 ft/sec 
e 

P2 = 50 psf 
H/6 = 0.283 



™total = E-4 slugs/sec 



AY 

E-4 

ft 



AH 

E-5 

ft 



M 






m 

E-5 

slugs/sec 



AH/AY 
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ALL ROWS ARE 
1/4" APART 
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The total area of the holes on one plate is larger than 
the calculated bleed area. This was done for three reasons: 

1) to bleed off a greater fraction of the boundary 
layer than calculated if needed (see figure (29)), 

2) to keep the same bleed area but vary the radial 
position of the effective bleed area (i.e., if 
bleed was desired at an inner radius, the outer 
holes were covered, or vice versa.), 

3) since the bleed plenum was fitted with pressure 
taps, flow velocities could be attained by 
securing the suction to the plenum and recording 
the pressure drop in the plenum due to the flow 
over the bleed holes. By covering all but one 
row of holes the velocity at a desired radius 
could be obtained by using the isentropic relation 
for pressure and Mach number or Bernoulli's 
equation . 

F. EXPERIMENTAL PROCEDURE 

Initially six nozzle rings were tested in the apparatus 
illustrated in figure (7) . Stagnation pressure applied to 
the nozzle stack, controlled by means of a regulator valve, 
was increased to a maximum of 130 pounds per square inch. 
Boundary layer bleed was attained by connecting the vacuum 
pump directly to the bleed plenums. The flow was observed 
using shadowgraph and Schlieren techniques. A typical run 
consisted of establishing the sensitivity of the Schlieren 
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(or shadowgraph) system to the desired level, starting the 
vacuum pump and allowing it to reach its lowest pressure, 
and then increasing the pressure to designated levels where 
photographs were taken. 

The number of nozzles was then changed from six to four 
and the spacing of the nozzle rings (the throat height) was 
increased from .05 inches to .09 inches. This change 
decreased the exit Mach number from 4.5 to 3.9. Since the 
total pressure loss across the shock was less, the back 
pressure did not need to be lowered to below atmospheric. 
Therefore the ejector nozzles were removed. Removing the 
ejectors allowed for longer runs at high pressures due to 
a decrease of approximately forty percent in the mass flow 
from the reservoir. 

The four nozzle rings were tested in the apparatus 
modified as shown in figure (32) in a fashion similar to that 
previously described for six nozzles. One important differ- 
ence was the use of a 200 cubic foot tank to increase the 
bleed capability of the vacuum pump. The pump was connected 
to the tank and was able to lower the pressure inside it to 
approximately one-half pound per square inch. The bleed 
plenums were then connected to the bleed tank and the suction 
was controlled with a gate valve. 

Flow conditions were observed and photographed at 
pressures up to the maximum attainable pressure of 250 pounds 
per square inch. Photographs were taken with and without 
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FIGURE (32) PHOTOGRAPH OF EXPERIMENTAL APPARATUS 



boundary layer bleed to compare the resulting flow fields. 
The bleed hole pattern was also altered to obtain pressures 
in the bleed plenum due to the velocity of the flow across 
the holes at various radii as discussed in Section III-E. 

A pointed probe was inserted into the flow at various 
locations and different total pressures to determine regions 
of subsonic and supersonic flow. 
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IV. RESULTS AND CONCLUSIONS 



A. RESULTS 

As was previously mentioned, six nozzle rings were 
initially installed in the test apparatus shown in figure 
(32) . Normal shocks were obtained in all except the end 
nozzles, but attempts at moving the shocks beyond the nozzle 
exit to form a "barrel" shock were unsuccessful at pressures 
of up to 130 psi.* Figures (33) (a) and (b) are Schlieren 
photographs of the six nozzles at pressures of 70 and 130 
psi, respectively. Before attempting higher pressures it 
was decided to increase the throat area in order to lower 
the exit Mach number and thus make it easier to move the 
shocks out of the nozzles. Since the sum of the throat areas 
must remain smaller than the feed area to prevent choking at 
the ends of the stack, the number of nozzle rings was reduced 
to four and the throat height increased to .09 inches. This 
resulted in an exit Mach number of 3.9 and increased the Mach 
number at the feed holes to .9. 

Three distinct flow patterns were observed in the flow 
over the range of pressures up to a maximum of 250 psi. 
Figures (34) , (35) and (38) are Schlieren photographs 

showing the development of the first stage flow structure. 
Figure (37) is a schematic of the flow just prior to its 



Applied pressures are gage pressures unless otherwise 
specified. 
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FIGURE (33) SCHLIEREN PHOTOGRAPHS OF NORMAL SHOCKS 
STANDING IN SIX NOZZLES 
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(a) NO FLOto 




FIGURE (34) SCHLIEREN PHOTOGRAPHS OF FOUR NOZZLES 
BEFORE "STARTING" 
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FIGURE (35) SCHLIEREN PHOTOGRAPHS OF CENTER THREE 
NOZZLES STARTING 
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(a) Flow at Pressure Less Than 50 PSI. 




(b) Resulting Coanda Jet Attached to Nozzle Wall 



FIGURE (36) SCHEMATIC OF COANDA JET FORMATION 
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FIGURE (37) SCHEMATIC OF FLOW AT A TOTAL PRESSURE 
OF 100 PSI. 
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FIGURE (38) SCHLIEREN PHOTOGRAPH OF FLOW JUST BEFORE 
TRANSITION TO SECOND FLOW STRUCTURE - 
TOTAL PRESSURE 125 PSI. 
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transition to the second stage at 145 psi. It can be 
readily observed that failure of the end nozzles (actually 
"half" nozzles) to start complicates the flow structure a 
great deal. Coanda jets appear to be attached to end nozzle 
walls creating a "bubble" of relatively slow flow. A slip 
layer results between this slow region and the Coanda jet 
and appears in the Schlieren photographs as a light area. 

It was originally considered that these bright regions could 
be expansion waves turning the flow back to each wall. 
Inserting a probe into this region as shown in figure (42-b) 
established that the flow was in fact subsonic since it 
failed to produce oblique shock waves from the tip of the 
probe. Pressure readings were also recorded on the wall at 
the nozzle exit and by assuming isentropic flow from there 
to ambient pressures (calculating upstream) , the Mach number 
was calculated and found to be subsonic. Table (5) lists 
the pressure data taken and the calculated Mach numbers. 
Figures (39) and (40) show the pressure and Mach number at 
the wall as functions of total pressure. The slip-line 
causes a bright area in the Schlieren photographs due to the 
higher entropy and temperature generated from friction, 
leading to an area of lower density. 

Careful analysis of the photographs of the flow up to 
pressures of approximately 145 psi. (e.g. before transition 
to the second phase) show a shock wave starting normal to 
the wall within the end nozzles and then curving to the 
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Table 5 Wall Static Pressure at Nozzle Exit and 
Calculated Mach Number. 



★ * 

Total Gage Wall Static Mach Number 

Pressure (psi) Pressure (psi) at Wall 



0 


12.94 






.43 






20 


12.42 






.50 






40 


11.59 






.59 






60 


11.0 






.66 






80 


10.51 






.71 






100 


9.72 






.79 






120 


9.28 






.84 






130 


- 


/ 


8.79 


- 


/ 


.89 


135 


- 


/ 


8.67 


- 


/ 


.90 


140 


8.56 


/ 


10.63 


.91 


/ 


.70 


145 


10.66 


/ 


10.73 


.69 


/ 


.69 


150 


10.69 


/ 


10.78 


.69 


/ 


.68 


155 


10.79 






.68 






160 


10.78 






.68 







The second values listed resulted from decreasing the total 
pressure from 150 psi (illustrating the hysteresis effect) . 
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TDTRL ERSE: PRESSURE CPSi) 

FIGURE (40) MACH NUMBER AT WALL VS. TOTAL PRESSURE 



angle of the nozzle wall, where the slip layer forms 



between the Coanda jet and the subsonic "bubble." The slip 
layer continues to curve back toward the end wall. Observable 
in the lower jet are two small, nearly normal shock waves — 
one where the shock attached to the wall begins to curve 
outward and one from the trailing edge of the nozzle ring. 
These combine to slow the flow to subsonic velocities as 
testing with the probe revealed. The jet on the upper 
nozzle is smaller and faster which is attributable to a 
slightly smaller opening at the throat. Figure (38) shows 
the alternately light and dark regions of a series of shock 
diamonds in the upper Coanda jet, indicating a higher Mach 
number than in the lower nozzle. The slip line in the upper 
case also bends back toward the wall. 

The high velocity jets occur for one of two reasons. 

First, since the end nozzles are really only "half-nozzles," 
the throats are only half the thickness of those in the 
center, in this case on the order of .045 inches. With 
these small openings located at the ends of the nozzle stack 
where the flow is entering axially at near sonic conditions, 
it has a difficult time "turning the corner" to exit out the 
first throat opening. If this were the case, whatever flow 
did make the 90° turn may attach itself to the nozzle wall 
as shown below. 
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A second, more feasible explanation is that regardless of 
whether or not the flow has trouble turning the corner, the 
flow that travels through the throat starts as a jet of air 
near the wall at low pressures (in the order of 50 psi.) is 
shown in figure (36-a) . This jet entrains the air near the 
nozzle wall and thus lowers the pressure in the wall region. 
The jet of air then moves toward the wall and attaches to 
it, displaying the "Coanda effect." As the flow in the jet 
moves out radially it slows and increases in pressure, 
forcing the slip line back toward the wall and creating a 
"bubble," observable in figure (35). As the total pressure 
increases toward the transition point of 145 psi. , the 
bubble becomes longer as the flow in the jet takes a longer 
distance to reach a higher pressure. 

The three nozzles in the center established normal 
shocks quite readily and moved outward to the nozzle exit 
with increasing pressure. As can be seen in the Schlieren 
photograph of figure (38) , connecting these three waves 
would result in a cylindrical shock bowed inward at the 
center. Reference 2 states that a loss of total pressure 
in an axial direction will manifest itself in a shock wave 
that is conical, tapering inward in the direction of 
decreasing total pressure, rather than strictly cylindrical. 
For a stack of nozzle rings fed from both ends, as was the 
case in this experiment, such a loss in total pressure would 
result in the condition observed at pressures between 90 and 
140 psi., that is, a tapering inward, from both ends. Total 
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pressures different from the region mentioned, however, did 
not result in a distortion of this type (see figures (41) 
through (48)). Seemingly, if axial total pressure losses 
were observable in a concave shock at pressures near 100 psi. 
they would also be observable at much higher pressures 
rather than disappear near 145 psi. 

At a pressure of approximately ten atmospheres, transition 
to the second phase flow structure occurred. Figures (41) (a) 
and (b) are Schlieren photographs of the flow at total 
pressures of 150 and 175 psi. This transition consisted of 

a sudden rise in the wall pressure, and therefore a drop in 

the Mach number, as shown in figures (39) and (40) , respec- 
tively. The Coanda jets no longer form a subsonic bubble, 
but instead the slip layer extends upward at an angle of 
19 degrees. The cause of the transition to the condition 
shown in figure (43) is the movement of the normal shock 
segments beyond the trailing edges of the nozzles. Since 
the flow leaves the nozzles at approximately a five to ten 
degree angle, oblique waves bend the ends of the normal 
shock segments outward giving them a concave shape. Shock 
diamonds form behind the oblique waves serving to return 

the flow to subsonic conditions. AVCO Everett in Ref. 2 

predicts that disturbances traveling downstream from the 
trailing edges would only affect the cylindrical shock wave 
locally, causing outward protrusions on the surface of the 
shock. Though the oblique shocks are rather large 
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FIGURE (41) SCHLIEREN PHOTOGRAPHS SHOWING CURVED SHOCK 
SEGMENTS BEYOND TRAILING EDGES OF NOZZLES 
AFTER TRANSITION TO SECOND PHASE 
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FIGURE (42) SCHLIEREN PHOTOGRAPHS WITH PROBE SHOWING 
REGIONS OF SUBSONIC FLOW 
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FIGURE (43) SCHEMATIC OF THE FLOW STRUCTURE AT A 
TOTAL PRESSURE OF 175 PSI. 
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FIGURE (44) SCHLIEREN PHOTOGRAPHS OF FLOW JUST BEFORE 
DEVELOPING INTO THIRD STAGE - TOTAL 
PRESSURE OF 200 PSI. 



92 




FIGURE (45) SCHLIEREN PHOTOGRAPHS SHOWING DEVELOPMENT 
OF FLOW TO THIRD STAGE 
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"disturbances," this prediction appears to hold by the fact 
that the ends of the normal shock segments "protrude" 
outward. Flow exiting in a more nearly horizontal direction 
would minimize the strength of these oblique waves, thereby 
making the shock segments less curved and closer to joining 
into a complete cylindrical shock wave. 

The second phase is also characterized by the disappear- 
ance of the subsonic "bubble" in each end nozzle. As the 
normal shock in the passage above (or below) the end nozzle 
moves beyond the trailing edge, a pressure mismatch occurs 
between the very low pressure region behind the shock and 
the Coanda jet coming off the trailing edge on the opposite 
side of the nozzle ring. This causes the Coanda jet to be 
diverted upward (or downward) , and the resulting slip-line 
to extend upward (or downward) . The angle that the slip-line 
forms increases from 19 degrees at transition (150 psi.) to 
29 degrees at the maximum pressure of 250 psi. This change 
is due to the increase in the pressure mismatch created at 
the trailing edge. Since the trailing edge is blunt, the 
Coanda effect will continue to force the jet to higher angles 
as it tries to move around the trailing edge. This is not 
the case with a sharp trailing edge where the angle of the 
Coanda jet would remain the same. 

As the pressure was increased from 150 psi. to 200 psi. 
the flow maintained the basic configuration shown in figure 
(43) . Figures (41) , (42) and (44) are Schlieren photographs 
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of the flow in this range of pressures. At 200 psi the 
oblique waves from the trailing edges are longer, better 
defined, and at a smaller angle. The curved shock segments 
have moved to approximately a quarter of an inch beyond the 
trailing edges. At approximately 210 psi. the flow structure 
is altered a second time. 

Transition to the third stage, unlike the abrupt change 
at 145 psi, , was a gradual development of normal shock segments 
in the nozzle wakes. Between pressures of 205 and 215 psi. 
the normal segments slowly appear, along with two weak oblique 
shocks at the trailing edges where one strong wave had been 
at lower pressures. Increasing the total pressure to 250 psi. 
changed the flow structure only by moving both the normal and 
curved segments of the shock slightly farther out from the 
trailing edges. Figure (47) is a schematic of the fully 
developed flow. Figure (46) is a Schlieren photograph of 
the flow at the maximum pressure applied, and shows that the 
normal portions of the shock are at a distance of approximately 
three eighths to a half of an inch. This distance is short 
of creating a lasing cavity wide enough to allow the theoret- 
ically available power in the flow to be extracted, which for 
four nozzles is 1.7 inches. For a spacing of the nozzles 
s , (s)/(r -r ) =0.85 , where r is the shock location. 

As a fraction of the minimum shock location allowable to 
extract the power available, (r -r ) , the location of the 
shock wave is: 
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FIGURE (46) SCHLIEREN PHOTOGRAPH OF FLOW AT MAXIMUM 
TOTAL PRESSURE OF 250 PSI. NORMAL SHOCK 
SEGMENTS APPROXIMATELY 3/8 INCH FROM 
TRAILING EDGES 



96 




FIGURE (47) SCHEMATIC OF THE FLOW STRUCTURE AT A 
TOTAL PRESSURE OF 250 PSI. 
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Boundary layer bleed was applied over the range of 
pressures with little if any effect on the flow. Figures 
(48) (a) and (b) are Schlieren photographs of the flow at 
200 and 250 psi. respectively with a one-half psi. bleed 
pressure applied to the first three rows of bleed holes on 
the upper plate. At a total pressure of 250 psi the flow 
entrained between the slip-lines is diverted toward the 
bleed plate. This was the only noticeable change caused 
by the bleed over the entire range of pressures. Attempts 
were made at applying enough suction to start the end 
nozzles but were unsuccessful. 

Homogeneity of the flow in the lasing cavity is a 
necessity for an operational laser. This was not the case 
for the fully developed flow at 250 psi. since the end 
nozzles failed to accelerate the flow to supersonic condi- 
tions. Other disturbances such as the oblique shock waves 
from the trailing edges should not be present to produce 
homogeneous conditions in the lasing cavity. 
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FIGURE (48) SCHLIEREN PHOTOGRAPHS SHOWING FLOW WITH 

BOUNDARY LAYER BLEED APPLIED TO UPPER PLATE 
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B. CONCLUSIONS AND RECOMMENDATIONS 

Testing of the cylindrical nozzle stack at a total 
pressure of 250 psi. established a shock wave approximately 
a half of an inch from the nozzle exit. The shock wave 
consisted of two segments normal to the flow (actually 
shock "bands" around the nozzle stack) connected by shock 
segments curved inward. Failure of the end "half-nozzles" 
to start resulted in subsonic flow near the end walls for 
the entire range of total pressures investigated. This 
prevented the normal shock from reaching the end walls. 

Since subsonic flow existed at the end walls and the 
shock did not reach the walls, the use of boundary layer 
bleed had very little visible effect on the flow structure. 

Flow visualization became very distorted at high 
pressures due to the extreme density gradients and the 
cylindrical geometry of the shock waves. 

The following recommendations are made for further 
research on this subject: 

1. Increase the diameter of the nozzle rings to allow 
for more mass flow to enter the nozzle stack without 
decreasing structural integrity and without choking 
in the feed passages. 

2. Decrease the thickness of the nozzle rings to lower 
the expansion ratio and thus lower the exit Mach 
number . 

3. Ensure sharp trailing edges on the nozzle rings. 
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4 . 



To help avoid starting problems in the end nozzles, 
use half of a nozzle ring (sliced along the edge) 
at each end. This would eliminate the "half-nozzle" 
arrangement and the Coanda jets. 

5. Investigate total pressures in excess of 250 psi. 

6. Install pressure taps on the center nozzle wall to 
determine Mach number. Optical distortion due to 
the large density gradients and the cylindrical 
geometry results in difficulty observing shock 
angles from a probe inserted in the flow. 
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APPENDIX A 



THEORETICAL LASING POWER AVAILABLE 

Intensity is defined as radiation power per unit area, 
with units of kilowatts per square foot. Assuming a 
Gaussian intensity distribution on each mirror as shown 
in figure (A-l) the intensity is 



I (r) = I 



-r 2 c 



To satisfy the boundary condition of 



I = I p (1/e 2 ) 



at r = r e , then 

o = 2/r e 2 . 



Therefore : 



I (r) 



I 



P 



-2r 2 /r 
' e 



2 



(A-l) 



Equation (A-l) yields a Gaussian peak centered at r=0. To 
obtain a ring with Gaussian cross section having the peak 
intensity at r=r, define 

r = r - r 
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FIGURE (A-l) LASING MIRROR INTENSITY LOADING 
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Inserting this into equation (A-l) yields 



- [E=£-yi] 2 

I(r) = I e r -r kwatts/cm 

To obtain the total power output, the intensity must be 
integrated over the mirror area. 

r m 2 

Power = f f I exp [- r r dr d0 

y~ v r 



P = 2 I / r exp[^— £-.y7 ] 2 dr . 

P r -r 

e e 



(A-2) 



An alternate method of finding the power available from the 
flow is 



Power = (mass flow rate) ( s P e cific^ 

energy 



Power = E (kjoules/lb^) m (lb m /sec) 



P = E m kwatts 



or 



P = 27ir G Lp 0 a 0 (p e /p 0 ) (a e /a c )M e E 



(A-3) 
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Equating equations (A— 2 ) and (A-3) : 



I / r exp [I-E-V2 - ] 2 dr = r L P a (p /p ) (a /a )M E , 



(A-4) 



where the only unknown in equation (A-4) is the outer radius 

of the mirror, r . 

' m 

In solving the integral Q on the left side of equation 
(A-4) define 



making 




/ 



and 



dx = 



V£_ 

r e~ F 



dr 




+ r . 



Making the corresponding changes in the limits of 
integration : 



@ r = r e , x = ~\J~2 

@ r = r m , x = -V2~ , 
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the integral now takes the form 



-VF 



r -r 2 r -r 

Q = / (r + x) e“ x ( ■■ - e ■•■ ) dx , 



VF 



VF 



VF 



or 



Q = 



-VF _ ¥ 2 _ F -VF 2 

/ r (!C) e" x dx + , / x e- x dx • 



VF 



VF 



VF YF 



(A-5) 



Since the second term on the right side of equation (A-5) 
contains an integral of an odd function between symmetric 
limits, the integral goes to zero leaving 



Q = 



VF 



■2 / 
0 



(r e -r) 

VF 



-x 



dx 



(A— 6 ) 



Multiplying and dividing the right side of equation (A-6) by 
it /2 and rearranging, 



Q = 



*V2tt‘ r (r-r ) 



, 2 -x , 

/ e dx 



o -yJV 



(A-7) 



Noting that the integral in equation (A-7) is the definition 
of the error function gives 



Q = V*/2 r (r-r ) erf(\fT ) 



(A-8) 
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Substituting equation (A-8) for the integral in equation 
(A-4) : 



I p [ \R/2 r (r-r e ) erf(\fT )] = V" p o a o {p e /p o } (a e /a o } M e E , 



or, solving for r , 



r = (r e /2) + 1/2 +C 



(A-9) 



where 



C = 



4 r L ME 

e e 

■\/tt/ 2 1^ erf (\[ 2 ~ ) 



Since 



r + r 

— e m 
r =5 / 



equation (A-9) can be written 



r m = r + C 
m e 



(A-10 ) 



The values in the expression for C are 



r =1.5 inches 
e 

p Q = (.002378) (32.174) (10) = .765 lb m /ft‘ 
a Q = 1119 ft/sec 
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.445 




(p / p o ) M=4.5 



.01745 



M 



e 



4.5 



E = 5 kwatts/ (lb^/sec) 

I =3 kwatts/cm^ = 2.8x10^ kwatts/ft^ 
P 

erf (/F ) = 2 ( . 9773 - .500) = .9546 . 



Using these values in the expression for C: 



C = .2692 (L) inches 



2 



For r g = 1.5 inches, equation (A-10) becomes 

r = V 2.25 + .2692 L , 

m r 

where L is the length of the nozzle stack in inches. 
Dimension L is illustrated in figure (7). 

As an example the power available from six nozzles was 
calculated. The theoretical lasing power available was 
found to be 



r 




r exp [ — ~\p2 } 2 dr 



r -r 
e 



= 2itI erf(\[2 ) \pi 7/2 [r(r-r )] 

P ^ 

= 2tt( 3) (30 . 5 2 /144) (.9546) (n/2) ** [1 . 635 ( . 135) ] 



Power = 32.16 kilowatts. 
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APPENDIX B 



CALCULATIONS FOR DETERMINING INLET (FEED) AREA 

Figure (13) shows the dimensions of the nozzle ring. 
The area of the hole segment between two centers is 

A = length x width 

= x (hole diameter) 

= 2tt. ( . 3/8 ^ 3 / - 4 i (1/8) x (3/8) 

A = .1407 in 2 . 

Since there are four segments: 



A 1 = 4 ( . 1407 ) 

A^ = .6627 in 2 . 



The area of the remaining half circles is 



a, = 8 x Liivsn, 



A 2 = .4418 in 



Therefore the total inlet area is 



A in “ A 1 + A 2 



A in = 1.104 in 2 . 
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APPENDIX C 



DETAILS OF THROAT HEIGHT CALCULATIONS 

The feed flow and the flow exiting through the nozzles 
must satifsy the continuity equation: 

6 

™in ~ ^nozzle 

n=i 

2 Pin u in A in “ 6 P* U * A * (C-l) 

where the factor 2 in equation (C-l) results from the 
nozzle stack being fed from both ends. 

Using the equation of state, and noting that velocity 
is the product of Mach number and the speed of sound, 
equation (C-l) can be written 



RT. ^ M in a in^ A in^ 3 (^r-) (1) (a*) (2nr t h) 



m 



P in P* 

— — M. A. = 6 — — h 
r~ m in r— 



Solving for the throat height: 



, _ , 1 , .. ,, , in. , T* . h 

h ^6 tt M in A in T . 

m 



(C-2 ) 
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P o T o t 

Multiplying the right side by (p— ) and (— ) 2 , 

o J ° 

equation (C-2) can be written as 



= (J-) 

'bT' 



o, /T*. % 



M in ( -F^ } ( P^ } ( t 
o o 



) 2 0 



o.h 

y • 

in 



(C-3 ) 



Using the isentropic equation relating stagnation temperature 
and pressure to static temperature and presure respectively, 
and applying them to the Mach equals one condition, 
equation (C-3) becomes 



h= 'Fir) M in '^X tftff)'- 8333 )''^ 2 -) 

o in 



or 



h = .0917 



M in (P i>/ p o> 



(C-4) 



Noting that 






equation (C-4) reduces to 



h = .0917 M. 

in 



( V T in>'‘* (T o /T in )Y/Y ' 1 
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or 



Y+l 

h = .0917 M. n (T in /T 0 ) Tt ^ rr) ' 



- (y+l) 

h = .0917 M in (l + I X" . 1 } .. M in 2 ) 2 (y-1 ^ 



h 



.0917 M. 
in 

( 1 + .2 M in 2 ) 3 
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APPENDIX D 



AXIAL VARIATION OF MACH NUMBER AND PRESSURE 



From Ref. 4, influence coefficients for flow in a 
channel are given as 




(1 



+ m 2 ) 



1-M 



dA 

A 



(D-l) 



and 



dP _ yM 2 dA 



(D-2) 



assuming isentropic flow and constant specific heat and 
molecular weight. Using the model shown in figure (D-l) 
and an inlet area of 1.105 square inches, 



A 



out 

A. 

in 



1 

IF 



out 



( rF } (A in ) 



Since 



A , = .06906 in 2 
out 



A 2 " A i - A out , 
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FIGURE (D-l) AXIAL FLOW MODEL 



A 2 ~ A 1 _ .06906 

h 1 1.10 5 - . 0697T6’ 



or 



dA 

A 



.06667 



Substituting this value into equations (D-l) and (D-2) and 
rearranging yields 



M i + i - M i - M i 



1 -M. 2 
1 



(.06667) ] 



and 



i+1 



Pi + P. 



Y ^ 2 

[ i-~ (.06667) ] ' . 

1 - M. 



(D-3) 



(D-4) 



A program was written for the Hewlett Packard model 9 830 -A 
computer to calculate the axial variation of Mach number and 
pressure using equations (D-3) and (D-4) . The program 
determines Mach number and pressure at stations one through 
thirteen using values of inlet Mach number and total 
pressure (PSF) as inputs. A listing of the program AXIAL 
and the results of the program are found in this appendix. 
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Table D-l 


Symbols for Program AXIAL 



Computer Notation 


Algebraic 


Description 


Ml 


M in 


Inlet Mach number 


P0 


P o 


Stagnation 
pressure (PSF) 


M( ) 




Mach number at 

subscripted 

station 


P( ) 




Static pressure 
at subscripted 
station 
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LISTING OF COMPUTER PROGRAM "AXIAL" 



:i P . i i l I; <Ji t Ni IT Mr! H# AND TO.P PRE . 

2: Dll'! i'P: ' 3 ] :■ Pi. 1 3 1 

3: INPUT N • F 

nm-Mi 

pij PC 1 ]”=P0 i. 1 +0 . 2*r1 1 •'!'£ I 5 > 

60 PRINT 

70 p R I i,!T " I HLET MACH#= " ; M 1 , " TOTAL PREP'S™ " 1 f 3 S P IP ' 
80 PRINT 



90 ! 


PRINT ‘ S'Tfi #" » " MACH# " 


t "P ST AT I C < PSP ) " 






100 


PRINT- 








1 1. 0 


PRINT " 1 " r. ME 1 3 5 PE 1 3 








120 


FOR I‘”2 TO 13 








130 


a™ i -nr i •••■]. ]t£ 








140 


ME I 3 -"ME 1-1 1 0.06667* 


ME 1-1 3* ( 1+0. 2* ME 


I 


- 1 3 !■? ) 


150 


PE I 3= PE 1-1 .] PPL 1-1 3*0 


. 06667* < 1 „ 4* ME I- 


; 


3 1 2 : ■ A 


160 


PRINT I > ME I 1? PE I 3 








170 


NEXT I 








180 


END 
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RESULTS OF PROGRAM 



"AXIAL" 



INLET MnCH#= 0 . i::. 



S'fH ¥f 


MRCTit: 


1 


6 - 2 




8. 13599 ■ 






0 „ 178065 .-70 




8. 1 6 1 89 9682 


J 


0. 1 508 1 5 754 


6 


9. 1 39737386 


r 


8. 130193914 




0. 1 21333930 


o, 


8. 1 13104282 


1 0 


0. 3.0544 6287 


1 1 


G. 698 521 528 


i 2 


0. 891689461 


1 3 


8.085514332 



INLET MflCH#= 0. 4 



Tfl # 


MAC 
0 . 


H# 

4 




0 „ 


36 629 6*5 7 


•j 


0. 


333172818 


4 


0. 


3 1 2 I 33254 


5 


0 , 


288627767 








b 


3 :i 


267286489 


i 


0. 


247321242 


C; 


8 . 


23080 1 659 


Q 


0 . 


21 3639648 


1 8 


0. 


198578371 


1 1 


0. 


184687414 


•| 


U . 


171852514 


1 


0 . 


1 59976780 



T ij ! H! ’-‘RE: 8 ! 

P STfiT I C PCr 

20378 . 0302 '? 

20658 ., 0539 

29727 ., 15727 ' 
20788 „ 891 83 
30933 , 5879 
28883 „ 38808 

20922 ., 18773 
28955 . 86257 
28985 - 891.84 
21818. 4726 
21832. 52282 
21851.68584 
21868-34491 



I 0 1 HI. P k E S 8 2 .!. 1 6 1 ' P I 

P ST AT I C < PSF > 

18951.28834 
19288. 12742 
1 9568 . 77096 
19884. 62462 
26884. 16139 
28 1. 7 3, 84 i 48 
20318. 7 1595 
28442 - 812 !. 1 
20549. 38963 
28641 = 1 1923 
28720. 21 889 
28738. 58753 
28347. 55658 
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RESULTS OF PROGRAM "AXIAL" (CON’T) 



H' T , 


r VCili-h 




TOTAL |:, PES 


TA # 




n fi c hi it 


P ST AT J C < P9F 






0 . 6 


1 . b . j y y u . j 4 -3 l ; j 


P 




G „ 5 3299 665 


1 7460 ,. 5 :!. 7 1. 4 






9 „ 400540045 


18107. 21786 


4 




0.1 ''"m-. ?o64 


18614. 67462 


.j 




L'l = 3995 7 7 0 : 11 ;) 


19624. 78217 


b 




0 „ b b o b 4 1 1 


1 9362 . 036 7 5 






O „ 3 7 t' 4 L" b • 0 


1 y b 4 3 ,. 1 4 b b j. 


o 




0 „ 1 1 0 0 4 ,: 'r G 3 


19879. 37194 


*;;! 




8. 288357130 


20S79. 237 lb 


i 0 




G« 2870481 50 


20249.28577 


i.t 




0. 247596104 


28394. 33385 


1 2 




8.229795217 


285 1 8 . 65806 


13 




8 . 2 1 3449337 


28625., 4283 



INLET MACHfe= 6=8 



TOTAL PRESS- 2 i. 166 



ST A # MriCH# P STATIC CPSF> 





8 . 8 


1 3 8 3 1 = 4 1 7 4 4 


y 


y . bo 28 b^bbb 


16184.31965 


o 


0 , bi c:: i f ,' r \ . 1 . , « , ! 


17194, 1841 


4 


9. 499706082 


1 7915= 50933 


5 


0,453085391. 


18472, 03545 


6 


0.413514623 


18917.40723 


•* 


8 , 379 1 2 1 795 


19231 , 61098 


8 


0, 348754325 


19583.71062 


9 


8. 321 639263 


19836, 82392 


18 


Q. 297226339 


28058 , 470 1 


11 


0.275105941 


20231,82379 


1 2 


8 . 254962478 


20386= 44623 


1 3 


8 . 2 -3 6 5 4 5 9 4 9 


28518=74135 
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APPENDIX E 



DETERMINING POINTS ON THE NOZZLE WAL L 

Following a development similar to Ref. 4, but changing 
the characteristic network as shown in figure (E-l) , 
equations can be developed to determine points along a 
right running characteristic proceeding in the -K direction. 
Along the characteristic segment 

3 a 3 

l ’Sri’ (v-0) = / sin y — ~ 6 dg (E-l) 

^ 2 

and 

3 a 3 

/ ( v— 0 ) d£ = - / sin y -- x g . 6 . d£ (E-2) 

where the negative sign on the right side of equation (E-2) 
results from traveling in the negative -£ direction when 
proceeding from point 1 to point 3. Assuming a very small 
characteristic mesh, the integrands on the right sides of 
equations (E-l) and (E-2) can be taken to be constant. 

With this assumption and some rearranging, the equations 
become 

v 3 = ? ^ v l +v 2^ + 

l sin 0 sin 0_ 

+ 7 [-sin y x — rj— A5 13 + sin y 2 — —1 An 23 ] 
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FIGURE (E-l) 



MODIFIED CHARACTERISTIC ELEMENT 



and 



e 3 = i ,v 1 -v 2 ) + £ (e 1+ e 2> 



, sin 0. sin 0 2 

+ j [-sin m 1 — AC 13 - sm y 2 — 



An 23 ] 
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Table E-l Results of Axisynunetric M. 0. C. 



MACH # 


V 


0 


V 


0+y 


0-y 


1.082 


67.57 


91 


1 




23.4 


1.256 


52.74 


95 


5 




42.3 


1.435 


44.17 


100 


10 




55.83 


1.605 


38.55 


105 


15 




66.45 


1.775 


34.29 


110 


20 




75.71 


1.95 


30.85 


115 


25 




84.153 


2.134 


27.945 


120 


30 




92.05 


2.14 


27.8 


122,65 


32.65 




94.85 


1.198 


56.6 


90 


3.55 


146.6 




1.35 


47.8 


93.57 


7.59 


141.37 


45.77 


1.53 


41.0 


48.36 


12.62 


139.36 


57.36 


1.70 


36.1 


103.23 


17.68 


139.33 


67.13 


1.873 


32.31 


108.12 


22.72 


140.43 


75.81 


2.05 


29.2 


113.03 


27.77 


142.23 


83.83 


2.24 


26.5 


117.95 


32.83 


144.45 


91.45 


2.349 


25.2 


120.58 


35.51 


145.78 


95.38 



1.49 


42.1 


90.0 


11.65 


132.1 


47.9 


1.66 


36.9 


94.75 


16.71 


131.65 


57.85 


1.835 


33.0 


99.58 


21.81 


132.58 


66.58 
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E-l Results of Axisyiranetric M. 0. C. (Continued) 



MACH # 


y 


e 


V 


6+y 


e-y 


2,02 


29.7 


114.43 


26.89 


134.13 


74.73 


2.21 


26.9 


109.30 


31.98 


136.20 


82.40 


2.415 


24.5 


114.18 


37.09 


138.68 


89.68 


2.52 


23.3 


116.78 


39.80 


140.08 


93.48 


1.835 


33.0 


90 


21.81 


123.0 




2.02 


29.65 


94.82 


26.95 


124.47 


65.17 


2.215 


26.85 


99.65 


32.08 


126.50 


72.80 


2.42 


24.40 


104.50 


37.22 


128.90 


80.1 


2.64 


22.20 


109.36 


42.40 


131.56 


87.16 


2.77 


21.15 


111.95 


45.15 


133.10 


90.80 


2.218 


26.75 


90 


32.18 


116.75 




2.42 


24.35 


94.81 


37.36 


119.16 


70.46 


2.65 


22.18 


99.64 


42.58 


121.82 


77.46 


2.90 


20.10 


104.48 


47.86 


124.58 


84.38 


3.045 


19.15 


107.06 


50.67 


126.21 


87.91 


2.655 


22.13 


90.0 


42.63 


112.13 




2.9 


20.1 


94.81 


47.95 


114.91 


74.71 
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E-l Results of Axisyminetric M. 0. C. (Continued) 



MACH # 


U 


0 


V 


6+y 


e-y 






" " ' ' 


' 


rr ' “ ' f 


3.19 


18.25 


99.63 


53.37 


117.88 


81.38 


3.36 


17.31 


102.20 


56.26 


119.51 


84.89 


3.19 


18.23 


90.00 


53.42 


108.23 




3.53 


16.45 


94.82 


59.03 


111.27 


78.37 


3.73 


15.55 


97.38 


62.03 


113.93 


81.83 


3.94 


14.70 


90.00 


64.98 


104.7 




4.19 


13.82 


92.55 


68.16 


106.37 


78.73 


4.47 


12.93 


90.00 


71.49 


102.93 




4.47 


12.93 


90.00 


71.49 


102.93 


77.07 


4.22 


13.70 


92.10 


68.66 


105.80 


78.40 


3.78 


15.35 


96.97 


62.68 


112.32 


81.62 


3.41 


16.92 


101.78 


57.01 


118.70 


84.86 


3.08 


18.95 


106.61 


51.41 


125.56 


87.66 


2.8 


20.9 


111.47 


45.87 


132.37 


90.57 


2.56 


23.0 


116.24 

119.92 


40.49 

36.24 


139.24 


93.24 
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Table E-l Results of Axisymraetric M. 0. C. (Continued) 



POINT 


MACH # 


y 


0 


74 


4.47 


12.94 


90 


63 


4.24 


13.65 


91 


60 


3.79 


15.32 


96 


49 


3.42 


17.03 


101 


46 


3.1 


18.85 


106 


35 


2.81 


20.8 


111 


32 






115 



76 


4.47 


12.93 


90 


65 


4.27 


13.54 


91 


62 


3.81 


15.22 


96 


51 


3.43 


16.96 


101 


48 


3.105 


18.79 


105 


37 


2.82 


20.8 


110 



78 


4.47 


12.93 


90 


67 


4.315 


13.41 


91 


64 


3.84 


15.11 


95 


53 


3.455 


16.8 


100 


50 


3.125 


18.63 


105 


80* 


3.5 




99 



* estimated 



V 


0+y 


6-y 


71.49 


102.93 


77.07 


68.87 


105.52 


78.22 


62.88 


112.04 


81.40 


57.18 


118.53 


84.47 


51.55 


125.16 


87.46 


45.97 


131.95 


90.35 


40.55 






71.49 


102.93 


77.07 


69.20 


105.05 


77.97 


63.15 


111.55 


81.11 


57.41 


118.04 


84.12 


51.73 


124.65 


87.07 


46.12 


131.47 


89.87 



71.49 






69.65 


104.44 


77.62 


63.59 


110.91 


80.69 


57.82 


117.29 


83.69 


52.10 


123.84 


86.58 



00 

87 

72 

5 

31 

15 

90 

00 

51 

33 

08 

86 

67 

00 

03 

80 

49 

21 
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APPENDIX F 



A. EJECTOR NOZZLE DESIGN 

To assure starting of the cylindrical nozzles a back 
pressure must exist which is lower than the total pressure 
downstream of the shock. For a stagnation pressure of ten 
atmospheres, the total pressure downstream of a normal 
shock at a Mach number of five is less than one atmosphere. 
Ejectors were designed to lower the back pressure to a value 
below the total pressure downstream of the shock assuming 
Mach-five flow. 

The models used for finding the flow conditions at the 
ejector nozzles are shown in figures (F-l) and (F-2) . 

Applying the continuity equation 

P 1 U 1 A 1 + p 2 U 2 A 2 = p 3 U 3 A 3 ' (F-l) 

the momentum equation 

P 1 U 1 2a 1 + P 1 A 1 + p 2 U 2 2a + P 2 A 2 = p 3 U 3 2a 3 + P 3 A 3 ' (F_2) 
and the energy equation 

P 1 U 1 A 1 ( S T 1 +U l 2/2) + p 2 U 2 A 2 (C p T 2 +U 2 2/2) = P 3 U 3 A 3 (C p T 3 +U 3 2/2) 

(F-3) 
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FIGURE (F-l) EJECTOR NOZZLE MODEL 
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D3 






FIGURE (F-2) TOP VIEW OF EJECTOR FLOW MODEL 
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to the control volume shown will .yield the flow conditions 
necessary at station one to assure the desired back 
pressure . 

The following assumptions are made: 



( 1 ) 


isentropic flow from station 4 


to 2 


( 2 ) 


isentropic flow from reservoir 


to ejector nozzle exit 


(3) 


P^ = 1 atm. 




(4) 


P T = 9.5 atm. 




(5) 


P = P 
*2 *1 




( 6 ) 


T Tl = 520 °R 




(7) 


M = 5.0 at shock 





Using normal shock relations the known flow at station (5) 
can be used to obtain conditions at station (4) . For 
M = 5.0, the pressure 

P 4 /P 5 =29.0 

P 4 = (29) (39.99) 

P 4 = .548 atm. , 

density 

P 4 /P 5 

P4 

p 4 
p 4 



= 5.0 

= (5.0) (p Q ) (P 5 /P 0 ) 

= (5.0) (.02378 g) (.01134) 

= .0434 lb /ft 3 , 
m ' 
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and temperature 



V T 5 



= 5.8 



T 4 = (5.8) (T q ) (T/T q ) 

= (518) (520) (.1667) 
T „ = 502.76 ° R 



can easily be found. 

Using the isentropic relations 

P 4 / p 2 = (P 2 /P 4 > = ( t 2 / T 4 ) ' 

all conditions at station two can be found in terras of the 
static pressure: 



p 2 p 4^ P 4 //P 2^ 



1/Y 



and 



p 2 = (.0434) (l|i£) 1/;L * 4 
P 2 



T 2 = T 4 (P 4 /P 2 ) (Y 1)/Y 



T 9 = (502,76) (-Jlii) 4 / 1 * 4 
2 P 2 



Since ™ 2 = ™ 4 , the velocity at 2 can also be found: 



m 



2 P 2 U 2^ 2A 2^ 



131 



where the area is taken as (2*A 2 ) since the control volume 
inlet area (A 2 ) is only one-half of the total exit area of 
the flow. Rearranging: 



m, 



u 2 " q 2 TT(a^-~a 1 )T * (144) 



_ (6.463) (144) 

" pJTT)T9'- .7S37 



U„ = 



55.9 



All conditions at station two are now known in terms of P 2 . 
Again using isentropic relations, conditions at station one 
can be determined in terms of : 



T-. -(y-D/y 

T = T (-pA) 

1 k' 1 



and 



T n = 



P-. = 



9 5 
520 (-£—■) 

1 



P 0 ^ ^ 



-.4/1.4 



’R 



9 5 - 1 / 1 - 4 3 

p ± = (.02378 g) (^-) lb m /ft 
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To find the velocity at station one: 



= T 1 (l + Iji M 1 2 ) . 



Since M, 



U 1 

yRTj^ ' 



u 



1 




T T 1 2 yRT, 

^ Y-l ^ g 



ft/sec 



Hence all conditions at station one are known in terms of 
P 1‘ 

If (=P 2 ) i s known, equations (E-l) and (E-2) can be 
combined to find and p^. A computer program was written 
to vary the values of D-^ and and to iteratively increase 
P^ from a set value until it satisfied the energy equation. 
An explanation of the program EJECTOR, a listing, and the 
results in tabular form are found in this appendix. Figure 
(26) shows the results in graphical form. The broken line 
at a pressure of 0.617 represents the total pressure down- 
stream of the shock assuming P Q (upstream) = 10 atm. and 
M = 5.0. Acceptable combinations of and must result 
in a pressure P 2 below this line. 
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B. EXHAUST LENGTH 



The length of the exhaust stack must be long enough to 
ensure complete mixing of the primary flow from the ejector 
nozzles and the secondary flow from the cylindrical nozzles. 

The model used for mixing of the two flows is shown 
in figure (F-3) . 

Bailey [Ref. 6] defines a non-dimensional spread 
parameter which is a good approximation to experimental data 



The non-dimensional similarity variable for the mixing 
zone was defined to be 




A value of ri = 1.84 gives the outer edge of the mixing zone 
where the velocity equals 99% of the centerline value. 

The minimum value of the mixing length is where the edge 
of the mixing zone touches the side walls, or when 




(F-4) 



where 



U 




mean 



7 



Y = Y 



max 
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FIGURE (F-3) EJECTOR FLOW MIXING MODEL 
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Substituting Y = Y and equation (F-4) into equation (F-5) 

msx 

the minimum mixing length. X . , can be found: 

mm 



mm 



crY 

max 

n 



= 24 



u 



mean 



U. 



< D 3 -°1 } 
2 (1.84) 



= 6.522 (4.85 - 1.75) 



(u 1 + u 2 ) 



2U. 



X m in “ ^ (1 + (3 - 1 » ' (F - 6) 

Using values from program EJECTOR for the combination of 
and chosen: 



U 1 = 1860 ft/sec 
U 2 = 664 ft/sec 



Putting these values in equation (F-6) : 



X min 



6.522 M . 664 . 

— (1 + IMTT ) 



(3.1) 



X ■ =13.72 inches 

mm 

The exhaust stack must therefore extend a minimum of 13.72 
inches beyond the ejector nozzle exit. Figure (F-4) shows 
the length of the exhaust stack with the ejector nozzles 
installed. 
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C. NOZZLE SHAPE 



For the chosen combination of = 1.75 inches and 
= 6.5 inches 

U x = 1860 ft/sec . 

(See section D of this Appendix.) This velocity is clearly 
supersonic. To determine the throat area required for the 
given exit diameter and flow velocity, the exit Mach number 
must first be found. Since 




M 



2 




/ 



and 



T = T q /(1 + 





/ 



then, after substituting and rearranging, the Mach number 
can be expressed as 



where 



M 




i 



K = 



u‘ 



yRgT 
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Inserting the appropriate values: 



M “ Vr^sr 

M = 2.49 . 



For isentropic flow through the nozzle 



A/ A* 



1 

M 



[{ y7T 



) (1 + 




M 2 ) ] 



Y+l 
2 (Y-l) 



For M = 2.49 



A/A* =2.61 . 



Solving for throat height: 



A* 



it (1 . 75) 2 , 1 , 

3 — ( 27ST ) 



A* = .921 in 2 , 



or 

Ua* 

D = v~ 

D* = 1.08 inches 

The ejector nozzles were fabricated from plexiglas with the 
throat diameter of 1.08 inches. Figure (F-4) shows the 
details of the ejector nozzle geometry. 



138 





FIGURE (F-4) EJECTOR NOZZLE GEOMETRY 
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D. EJECTOR 



The program EJECTOR solves the continuity and momentum 
equations for different combinations of nozzle diameter and 
exhaust stack width assuming a starting value of pressure 
at stations one and two of 0.37 atmospheres. The pressure 
is increased in increments of 0.0001 atmospheres until 
the flow conditions at 1, 2 and 3 satisfy the energy equation. 
When the energy equation is satisfied, the values are printed, 
and the process is repeated with new values of D-^ and . 
Inputs to the program include conditions immediately 

3 

downstream of the shock: (lb m /ft ) , T^ (°R) , 

(atm.) and m^ (lb^/ft ) ; and total temperature (°R) 

and total pressure (atm.) at station one. 

Outputs are and (both in inches) , P-^ (=P ^ ) in 
atmospheres, and and , both in ft/sec. 
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Table F-l Symbols for Program EJECTOR 



Computer Notation Algebraic 



A A2 A3 
PI P2 P3 P4 
R1 R2 R3 R4 
Ml M2 M3 M4 
G 
R 

Cl 



A 1 A 2 A 3 




Description 

2 

Flow area (in ) 

Static pressure (atm.) 

Static density (lb m /ft^) 

Mass flow (lb /sec) 
nv 

Specific heat ratio 

Gas constant for air 
(53.3 ft-lb f /lb m °R) 

Specific heat @ constant 
pressure 



J 

Tl T2 T3 T4 

T 

P 

U1 U2 U3 

K 

L 

Dl 

D3 



T 1 T 2 T 3 T 4 



Tl 



Ti 



U 1 U 2 U 3 



778.26 ft-lb f /BTU 

Static temperatures (°R) 

Stagnation temperature (°R) 

Stagnation pressure (atm.) 

Velocities at stations 1, 2 
and 3 respectively (ft/sec) 



2*g 



C *J 
P 

Diameter of one of the two 
ejector nozzles (inches) 



D 



3 



Length of one side of exhaust 
stack . 
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D. LIST OF COMPUTER PROGRAM "EJECTOR" AND RESULTS 



1 : : N T ,NPU i -Or!U I'i , ON’: H STB. 

!. : I r * UpU T PHO'KLE; *U F‘ KHK'i'. 





. 1. 1 i ! ' U i r 


• *•'■ '* ! h i" 4 ! i M'r 






- : 


4 : 


ji'iPlP! ! 3 T ! 


i! K 


414 t 


1 . ONs •• 4' . 1 IT : 




PRIM! 


INPUT' "f T 1 ( 


. R > , P 


'TH 


Pi I Pi ) " 


b 0 


i s PUT ■ 


? F 








I d5 


i 










120 


G- 1 . 4 










. 1 . 30 


! • j J .. : 










140 


• «... -? 

•J f I * 


8 C 








.150 


C ! ~ 8 „ 


.4 I 








I so 


PR • N T 


"IK " « .03" :■ 


" P 1 " 


U 


! 1 ' " UP " 


1 0 


i ■' ft 1 i' i i 










! 74 


FOR H 


-1 TO 2.5 


STEP 


0 „ 




1 76 


ft =2 *3. 


1 4ibau 1 2.- 


'(4*1 


44 ) 




i 77 


P 1 =0 „ 3 










1 78 


FOR 171 


-6 TO 9 STEP 8 


. 5 




1 88 


HOLDS'! 


2/ 144 








1 8 b 


H2 = H3 


hi " • '• 4! * ( 3 •• ' J 


>*i)i 


••••' 1 44-Pi 2 > 


200 


P2 = P 1 










2 1 0 


R2=R4 > 


( P4/R2 > 1 1 i 


/G > 






228 


T £'•="!' 4 


<P4.-"P£)ta 


G- i ) 


/ G ) 




230 


U2=M4/ 


( R2*fi2 0 








2 4 8 


Tl-T*< 


R/p i > t < i 


G ) G 


'i 




208 


P 1 = ( P-" 


1 P* i ) ) ■ ■ ( r •••' 


P 1 > t 


< i 


•■4.:, 1 b 


268 


U1=SQR 


( ( S . T ] > 1 


> *2* 


Li & K 


*11*32. 17 ■ 


278 


i33=R 1 * 


U 1 * Pi H 3 + K 2 


Hi 2* 


82/ 


83 


2 88 


U3= ( R 1 


4Ul'T' t :.. i-H/liS 


•HP 2 


....pc, 


1*21 lb* 22 1 7 + P,:i-' 


290 


R3~P 1 * 


U 1 *H- / (. U3*PI 


3 > -t R 


2*U 


2 * r! 2 ( U 3 * H 3 1 


3 0 8 


T3-P3* 


2 1 2'. c- •••'' ( K 4 P 3 








3 1 8 


M 2 ~R 1 * 


U .1 4 PI 








328 


?-i2-R2:- 


i.l 2 4 Pi 2 








330 


M3-R3* 


U3*R3 








348 


!<=2*32 


. 1 7 








35 8 


L- : C 1 * J 










368 


Hi =M1* 


( L. * T 2 -i i 2 1 1 2 


••••' K ) 






3 78 


U2--M 


( L*7 2+U2'T2 


/!< > 






3 8 8 


Pi Hi 1 +U 










390 


X=M3* < 


L*T'3fU3-t'2.-' 


;• 






400 


d=w -x 










482 


IF IK 1 


O THEN 410 








483 


if p i :> 


0.65 THEN 


4 i 8 






404 


p 1 =p 1 + 


0 . 80 1 








4 0 6 


GOTO 2 


0 0 








4- :i i:3 


PRINT 


U 1 ... 113 J P 1 .. I.J 


j. 5 i J2 






428 


HFXT Ii 


3 








438 


PRINT 










4 48 


NEXT I) 


1 








458 


END 
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